Exploring the anomaly in the interaction cross section and matter radius of ^^O 
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New measurements of the interaction cross sections of ^^-^^q goOA MeV performed at the 
GSI, Darmstadt are reported that address the unsolved puzzle of the large cross section previously 
observed for ^"^O. The matter radii for these oxygen isotopes extracted through a Glauber model 
analysis are in good agreement with the new predictions of the ab initio coupled-cluster theory 
reported here. They are consistent with a ■^■^O+neutron description of ^'^O as well. 



PACS numbers: 25.60.-t, 25.60.Dz, 21.10.Gv, 21.60.De 

The nucleoli magic numbers are the fundamental ba- 
sis for the concept of nucleons being arranged in a shell 
structure. While the distributions of nucleons in stable 
nuclei are fairly well understood, the neutron-rich nuclei 
show signatures of unconventional behaviour. Of partic- 
ular interest is the region around the new magic number 
iV=16 [IHS], at the neutron drip-line that is a new bench- 
mark point. Here, the nucleus ^'^O plays a special role 
as a large enhancement in the interaction cross section 
was observed for ^^N, ^■^O and ^^F [1 . This leads to 
a large matter radius which was thought to reflect the 
formation of a neutron halo. However, the prediction 
by a core(^^0)-l-(2si/2) neutron halo model was much 
lower than the data, and the matter distribution of ^^O 
remained an unsolved puzzle being a challenge for all 
models of nuclear structure. 

The purpose of this Letter is to address this crucial 
issue, by measuring the interaction cross section and re- 
liably extracting the root mean square radii, in order to 
reach a conclusive understanding on ^^O. The new obser- 
vations resolve the existing anomaly, showing a smaller 
cross section consistent with both new ah initio model 
predictions presented here, and a ^^O core-|- neutron de- 



scription. 

There are several indications of a sub-shell gap at 
iV=14 in the oxygen isotopes. First, the 2+ excitation 
energy of 3.199(8) MeV in "^^O [1 [5] is large. Second, 
proton inelastic scattering of ^^O [B] revealed a small de- 
formation parameter /3=0.26±0.04. This sub-shell gap 
suggests that the ^'^O ground state could have a large 
component of single-particle configuration with 2si/2 va- 
lence neutron. The nuclear and Coulomb breakup mea- 
surements of one-neutron removal from ^^O reported 
the 2si/2 spectroscopic factors of 0.97±0.19 [3 [5] and 
0.78±0.13 |9j respectively. Significant yield of ^^O ex- 
cited state (3.2, 5.8 MeV) components were observed 
in the nuclear breakup. However, the ^^0-|-n(2si/2) de- 
scription in which the ^^O core is considered identical to 
the bare ^^O nucleus severely underpredicted the mea- 
sured interaction cross section T]. It was therefore pro- 
posed [10^ that the ^^O core within ^"^O is possibly mod- 
ified and enlarged compared to the bare ^^O nucleus, 
giving rise to a larger interaction cross section. The rela- 
tively narrow momentum distribution of the two-neutron 
removal fragment ^^O from ^^O pJJ suggested that ^^O 
might have some probability of two neutrons occupying 
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the 2si/2 orbital. However, this is not consistent with 
the reported 2si/2 spectroscopic factors. 

The energy gap at N—IA in ^-^O was found to be 
2.79(13) MeV from fragmentation of ^^Ne, populating the 
resonance m at 45(2) keV above the 220+n thresh- 
old [12]. This was considered to be the 5/2+ excited state 
in ^'^O. The higher lying resonance at 1.3 MeV above the 
neutron threshold observed in the ^^0{d,p) reaction [T3] 
was understood as the 3/2+ excited state that shows a 
4.00(2) MeV gap between the 2si/2 and 1^3/2 orbitals. 

Recently, the neutron knockout of ^''O — ?>^'^0 revealed 
a large two-neutron spectroscopic factor of 1.74±0.19 ex- 
hausting the s-orbital occupancy with no significant ob- 
servable d- component establishing a spherical shell clo- 
sure at iV=16 [Umi. The neutron-unbound excited state 
of at 4.72(11) MeV [31, considered to be the first 2+ 
state, is also supportive of ^^O being a doubly-magic nu- 
cleus. These findings are difhcult to reconcile with the 
unusually large interaction cross section of ^'^O reported 
in [T]. Such an anomaly may point to unknown struc- 
ture effects to which the breakup reactions may not be 
sensitive. The issue therefore is of utmost importance to 
be addressed since the neutron-rich oxygen isotopes form 
crucial benchmark points for understanding the evolution 
of shell structure in neutron-rich regions. It has been 
shown that three-body forces play an important role to 
define the drip-line of the oxygen isotopes [TS] . 

We performed an experiment to measure the inter- 
action cross sections of 22,23q fragment sep- 
arator FRS at GSI [T^. The experiment layout is 
shown in Fig.l of Ref.[T3- The measurement is done 
by the method of transmission where the total inter- 
action cross section for reactions in the target is given 
by cr/=(— l/i)ln(i?i„/i?o„t). The transmission ratio is 
Rin = n/^/N-^ where N-^^ and Nj^ are the number of 
"^O before and after the target, respectively. Rout is the 
same but for an empty target and t is the number of 
target nuclei per unit area. 

The '^O nuclei were produced from the fragmentation 
of a lA GeV '''^Ca beam interacting with a 6.3 g/cm^ 
thick Be target. The fragmentation products were sepa- 
rated and identified using the first half of the FRS, where 
plastic scintillator detectors placed at the two dispersive 
foci, Fl and F2, measured the time-of-flight (TOF). Pre- 
cise beam position measured using time projection cham- 
bers (TPC) detectors placed at F2 and the magnetic 
rigidity of the FRS provide information on the mass to 
charge ratio {A/Q) of the fragments. A multi-sampling 
ionization chamber (MUSIC) placed at F2 measured the 
energy-loss (with a la resolution of ^ 3%) providing the 
Z identification of the '^O beam. 

A 4.046 g/cm^ thick carbon reaction target was placed 
at F2. The second half of the fragment separator consists 
of a dispersive focus at F3 and an achromatic focus at F4. 
The ion optics mode was selected to match the dispersion 
of the first half with the second half. The magnetic rigid- 
ity of the second half was set to transport the unreacted 
to the final focus (F4). Here products were identi- 



1400 




900 [ I I I I I I I I I 

15 16 17 18 19 20 21 22 23 24 25 

Mass Number (A) 

FIG. 1: The interaction cross section of '^O+C as a function 
of the mass number. The circles represent the present data, 
and the squares are from Ref. 'T . The line shows the A^^^ 
dependence normalized to ^"O. 



tied in A/Q using magnetic rigidity, TOF between plastic 
scintillators at F2 and F4, and position measurement us- 
ing TPC detectors placed at F4. Two MUSIC detectors 
were placed at F4 to measure the energy-loss of the prod- 
ucts after the reaction target. Events that were observed 
to be consistent with Z—8 in either MUSIC were counted 
as unreacted events. To account for losses occurring due 
to reactions in other materials in the setup, data were 
also collected with an empty target frame. The position 
and angle of the incident beam before the target at F2 
were restricted from beam tracking measurements such 
that transmission from F2 to F4 was a constant within 
the selected phase space. 

The measured interaction cross section at ^900A MeV 
is shown in Fig. 1 as a function of the mass number of 
the oxygen isotopes with the filled circles showing the 
present data from Table 1. The open squares are the 
data of Ref.[T]. It is seen that the present data of ^'^O is 
smaller than that reported earlier. The total uncertain- 
ties shown are dominated by statistics. The uncertainty 
from contamination of the incoming beam (for Z and A) 
was at most ~ 0.9%, and that from the target thickness 
was 1.2 %. Transmission uncertainty was ~ 1.5%. The 
cross section for ^^O is also slightly smaller than previ- 
ously reported. The cross section for the isotope ^^Ne ex- 
tracted from the same data results in 1271±35 mb which 
agrees within uncertainty with the value of 1244±44 mb 
reported in Ref. [18 . The interaction cross section of ^■^O 
reported here is ~ 8-9% larger than ^^O, which may not 
be sufficient to be classified as a one-neutron halo (^^Be, 
a halo, has ~ 16% larger cross-section than ^'^Be). This is 
consistent with the fact that the one-neutron separation 
energy of ^^O, 5'„=2.7 MeV is quite large which inhibits 
the tunnelling of the wavefunction into the classically for- 
bidden region to form a halo. 

The root mean square (r.m.s.) matter radius is ex- 
tracted by interpreting the data in the framework of 
the Glauber model including the higher order terms that 
are missing in the usual optical approximation [191 120] . 
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The matter density was considered to be a Fermi den- 
sity of the form p{r) = pq/{1 + exp((r — R)/a)) where 
R — rf^A^/^. The calculated cross sections for differ- 
ent values of radius (ro ~ 0.8 — 1.4 fm ) and diffuseness 
(a ~ 0.3 — 0.7 fm) parameters are shown by the different 
filled points in Figs. 2a and 2b. Calculations using sepa- 
rate neutron and proton Fermi functions were found to 
yield similar r.m.s. radii. The r.m.s. matter radius that 
can reproduce the experimental cross section is found to 
be 2.75±0.15 fm for "^^O while for a radius of 2.95± 
0.23 fm is extracted (Table 1, shown by the horizontal 
arrows in Figs. 2a and b). To investigate any model de- 
pendence we also extract the r.m.s. radii with a harmonic 
oscillator form of density [18] . The radii extracted shown 
by solid line (open points) in Figs. 2a and 2b are con- 
sistent with those using Fermi density (Table 1). Since 
the oscillator width is the only one parameter here, the 
uncertainty of the radius is found to be slightly smaller. 

A sub-shell closure at A^=14 has been discussed for ^^O 
which allows us in a simplistic model to describe ^^O be- 
ing composed of a ^^O core-f neutron and we interpret 
next the interaction cross section in a few-body Glauber 
model framework [21]. The core ^^O is considered to 
have the same Fermi density profile as the bare ^^O nu- 
cleus mentioned above. The wavefunctions of the valence 
neutron are calculated with a Woods-Saxon bound state 
potential (ro=1.27 fm and a=0.67 fm) where the depth 
is varied to reproduce the effective neutron separation 
energy. The dashed line in Fig. 2c shows the calculated 
cross sections with the valence neutron in the 2si/2 or- 
bital. The cross section shown is calculated as a function 
of the rms matter radius of ^^O (from Table 1). The 
horizontal shaded area is the measured interaction cross 
section of ■^'^O. The overlap of the calculated values with 
this shaded region shows consistency with data. It is seen 
therefore, that ^■^O can be described by a ^^O core-|-2si/2 
neutron. The dotted line shows a similar calculation but 
for the neutron in the ld5/2 orbital and with the core 
^^O in its 2+ excited state. The density of ^^O in its 2+ 
state was assumed to be the same as the ground state. 

To gain a better understanding, we perform ab initio 
coupled-cluster (CC) [12] computations. This method 
is uniquely suited for describing nuclei with closed neu- 
tron and proton subshells and their neighbors. Within 
particle-removed CC theory, the ground state of ^"^O is 
described as a superposition of Ih and lp2h excitations 
on top of the correlated ground state of ^^O We 
employ a low-momentum version of the nucleon-nucleon 
interaction from chiral effective field theory [H] that 
results from a similarity renormalization group transfor- 
mation |25j and is characterized by a momentum cutoff 
A. We work with the intrinsic Hamiltonian, i.e. the ki- 
netic energy of the center of mass is subtracted from the 
total kinetic energy. As a result, the coupled-cluster wave 
function factorizes with a high degree of accuracy into a 
product of an intrinsic wave function and a Gaussian for 
the center of mass 26J. In this framework we compute 
the density of the ^'^O ground state. The intrinsic den- 
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FIG. 2: The interaction cross section data for (a) ^^0-|-C and 
(b) ^^0-|-C (shaded regions). The squares are Glauber model 
calculations with different Fermi density parameters ro and 
a shown as a function of RI,'J„s('''0). The dotted lines show a 
guide to the eye for determining the limits of consistent 
with the data. The solid line are calculations using harmonic 
oscillator density, (c) The interaction cross section data for 
^^O-l-C (shaded region). The dashed (dotted) line is ^^O+n 
few-body Glauber model calculation for 2si/2 (ld5/2) orbitals 
as a function of different r.m.s. radii within uncertainty 
of the value quoted in Table 1. 



sity, i.e. the density with respect to the center of mass, 
results from a deconvolution with respect to the Gaus- 
sian center-of-mass wave function. It enters the com- 
putation of the interaction cross section and the matter 
radii. Our model space consisted of 30 bound and con- 
tinuum Woods-Saxon orbitals for the neutron, 1=0 and 
1=2 partial waves, and 17 major oscillator shells for the 
remaining neutron partial waves and the protons. 

The computed point matter radii are in good agree- 
ment with our data (Fig. 3a). Smaller values of the mo- 
mentum cutoff lead to smaller radii. The relative uncer- 
tainties in the experimental radii are larger than those 
of the interaction cross sections due to the uncertainties 
in the Fermi density profiles. We also compare the cross 
sections calculated with density distributions from the 
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FIG. 3: (a) The red filled circles show R™^'^'^'"''' , dashed ver- 
tical line is R™^"' and the blue open squares are R™^ from 
[1]. The diamonds /squares / triangles are coupled-cluster 
calculations with cut-off parameter = 4.0 / 3.8 / 3.6 fm~^ 
where the filled points are R™s Eind the open ones are pro- 
ton r.m.s. radii respectively, (b) The interaction cross section 
of -^O-fC as a function of the mass number. The red filled 
circles are present data, the squares are data from Ref. '1 . 
The diamonds /squares / triangles are CC calculations with 
cut-off parameter — 4.0 / 3.8 / 3.6 fm"'^. The cross marks 
(open circles) show results from Ref. [27] (Ref. [28]). 

coupled cluster theory in Fig. 3b. The cross sections with 
A — 4.0 & 3.8 fm"^ are in good agreement with the data, 
while that with A=3.6 fm"^ is slightly below the la er- 
ror. The variation of the radii with the cutoff provides an 
estimate of the contributions of neglected short-ranged 
three-nucleon forces. The relative isotopic differences in 
radii from ^^~^*0, depend only very weakly on the cut- 
off, suggesting it to probably have a weak dependence 
on three-nucleon forces. The results with densities from 
Ref. |27| which were generated by a Slater determinant in 
a mean field potential are shown by the cross marks in 
(Fig. 3). The open circles (Fig. 3) show results with 
densities from a Skyrme-Hartree-Fock potential [28] . 
In conclusion, this work reports new measurements of 



the interaction cross sections of 22,23q gQQA MeV. The 
new data for ^'^O is consistent, within experimental er- 
rors, with a model of ^^O core -I- valence neutron in the 
2si/2 orbital thereby addressing the existing anomaly in 
its structure. Beyond this simplistic model, new coupled- 
cluster calculations reported here, show excellent agree- 
ment with the present data. This shows the significant 
advancement jointly by experiment and ab initio theo- 
ries in reaching a conclusive understanding on the mat- 
ter distribution in ^'^O. The radii extracted are shown to 
be consistent for the different parametric density forms. 
The results show growth of neutron skin from ^^O to ^"^O, 
but a large enhancement characteristic of a halo is not 
observed for ^■^O. 
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Isotope 


ai(Acr) Acr(Stat.) Acr(Syst.) iJ^if"™'' 
(mb) (mb) (mb) (fm) 


^rms 

(fm) 


230 
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1216(41) 33.1 24.7 2.95(0.23) 


2.75(0.07) 
2.97(0.11) 
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